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Chapter  I.  Introduction 


Electromagnetic  wave  propagation  in  conducting  media  has  been  of  practical 
interest  since  the  beginning  of  the  century.  Toward  the  end  of  the  first  World  War, 
limited  experimental  and  theoretical  work  focused  on  the  generation  of  electromag¬ 
netic  fields  in,  on,  and  above  the  sea  by  submerged  cables  carrying  alternating  cur¬ 
rent  [Drysdale,  1924;  Butterworth,  1924].  This  work  was  motivated  by  the  use  of 
cable-generated  electromagnetic  fields  for  navigation  [Wright,  1953].  These  fields  are 
also  of  interest  because  of  the  possibility  that  they  could  be  used  for  communicat¬ 
ing  with  submarines  [Moore,  1951,  1967].  Despite  the  absorption  of  electromag¬ 
netic  energy  in  sea  water,  electromagnetic  signals  with  frequencies  in  the  ULF /ELF 
range  (freqencies  less  than  3  kHz)  are  able  to  propagate  to  moderately  great  depths 
in  the  ocean  and,  as  a  result,  they  can  provide  a  means  of  communication  with 
deeply  submerged  submarines  [Wait,  1972].  Another  application  of  these  signals  is  in 
geophysical  prospecting.  A  significant  amount  of  research  has  used  electromagnetic 
techniques  to  study  the  structure  of  the  Earth’s  crust  [Burrows,  1963;  Wait  and 
Spies,  1972a, b,c].  Propagation  of  ULF/ELF  signals  through  the  earth  has  also  been 
of  considerable  interest  for  mine  communication;  in  the  event  of  a  mine  disaster, 
telephone  wires  and  other  normal  links  of  communication  could  be  interrupted,  but 
communication  with  the  trapped  miners  could  still  be  achieved  through  the  Earth 
[Wait  and  Spies,  1973]. 

The  first  extensive  theoretical  study  of  electromagnetic  wave  propagation  be¬ 
tween  submerged  stations  in  sea  water  was  reported  by  Moore  [1951],  and  most  of 
the  subsequent  theoretical  work  centered  on  a  sea  of  infinite  extent  [Wait,  1952a; 
Kraichman,  1970]  or  on  the  propagation  effects  associated  with  the  sea  surface 
above  a  very  deep  sea  [Moore,  1951;  Wait  and  Campbell,  1953;  Hansen,  1963; 
Banos,  1966;  Fraser-Smith  and  Bubenik,  1976;  Bannister  and  Dube,  1977],  These 
sea-surface  studies  indicated  that  if  the  horizontal  distance  between  the  source 
and  observer  is  larger  than  several  skin  depths,  the  energy  received  may  follow  an 
up-over-down  mode  above  the  surface  [Moore  and  Blair,  1961].  This  result  led  to 
the  conclusion  that  near  the  sea/sea-bed  interface,  the  energy  received  may  follow 
an  analogous  down-under-up  mode  in  a  weakly  conducting  sea  bed  under  certain 
conditions  [Bubenik  and  Fraser-Smith,  1978]. 


Such  a  mode  has  been  investigated  in  a  somewhat  different  context  by  a  number 
of  researchers  [Wheeler,  I960;  Burrows,  1963;  Mott  and  Biggs,  1963;  Wait  and  Spies, 
1972a.b  c  King  and  Smith,  1981],  These  studies  suggested  that  a  waveguide  may 
exist  under  both  the  sea  and  the  continents  and  that  it  may  become  a  usable 
communication  link  if  other  above  ground  communications  are  disrupted;  it  would 
also  have  much  lower  noise  levels.  In  an  electrically  shallow  sea,  the  sea-surface 
and  sea-floor  modes  may  exist  simultaneously.  Weaver  [1967]  analyzed  the  fields 
of  electric  dipoles  submerged  in  sea  water,  and  Ramaswamy  [1972]  considered  a 
submerged  horizontal  magnetic  dipole,  both  taking  a  sea  that  was  one  skin  depth 
deep.  Ramaswamy  also  computed  the  field  components  when  both  the  source  and 
receiver  are  located  at  the  sea/sea-bed  interface  as  a  function  of  sea  water  ind  uction 
number  varying  from  0  to  3  Numerical  data  describing  the  electromagnetic  fields 
produced  by  a  vertical  magnetic  dipole  submerged  in  a  sea  of  finite  depth  [Bubenik 
and  Fraser-Smith,  1978;  Fraser-Smith  and  Bubenik,  1979)  reveal  the  effects  of  a 
strongly  conducting  sea  bed  Other  electromagnetic  methods  for  obtaining  sea¬ 
bed  conductivities  have  been  proposed  [Brock-Nannestad,  1965;  Bannister,  1968a, b; 
Coggon  and  Morrison,  1970).  Bannister  [1968 a]  observed  that  the  conductivity  of 
the  sea  bed  may  be  determined  by  measuring  only  the  horizontal  component  of  the 
magnetic  field  produced  at  the  sea/sea-bed  interface  by  a  long  horizontal  line  source 
located  at  the  air/sea  interface  Coggon  and  Morrison  [1970]  considered  a  vertical 
magnetic  dipole  submerged  in  sea  water  and  analyzed  the  electromagnetic  fields 
over  various  ranges  of  sea-water  induction  numbers  and  sea-bed  conductivities.  A 
recent  active-source  electromagnetic  sounding  experiment  on  the  ocean  floor  [Young 
and  Cox,  1981]  introduced  new  possibilities  for  studying  the  electrical  conductivity 
of  the  ocean  crust. 

This  work  considers  an  electrically  deep  sea  and  calculates  the  electromagnetic 
fields  produced  on  the  sea  floor  by  three  different  sources  co-located,  along  with  the 
receiver,  on  the  sea  floor,  as  illustrated  in  Fig.  1.1.  The  effects  of  the  sea  surface  are 
neglected  because  attenuation  along  the  paths  of  propagation  to  and  from  the  sea 
surface  is  large.  Both  the  sea  and  sea  bed  are  assumed  to  be  isotropic,  homogeneous, 
and  time-invariant  conducting  media  separated  by  a  plane  interface.  Computations 
are  restricted  to  frequencies  in  the  ULF/ELF  range,  since  it  is  only  signals  in  these 


bands  that  can  propagate  to  significant  distances  through  the  conducting  media 
before  becoming  severely  attenuated.  The  displacement  current  terms  in  both  media 
are  neglected,  and  this  assumption  is  well  justified  in  these  frequency  ranges.  It  is 
assumed  that  both  media  are  nonmagnetic,  with  permeabilities  equal  to  the  free* 
space  permeability  (/io  =  4jt  X  10-7H/m).  The  conductivity  of  the  sea  water  is 
it,  =  4  S/m  and  the  conductivity  of  the  sea  bed  varies. 

The  three  sources  considered  are  the  vertical  magnetic  dipole  (VMD),  horizontal 
electric  dipole  (HED),  and  infinite  cable  (IC).  It  is  assumed  that  all  the  fields  vary 
with  time  as  exp(iu?0-  The  amplitudes  of  the  field  components  for  the  VMD  and 
HED  are  obtained  numerically  through  the  techniques  described  by  Bubenik  [1977], 
Three  of  the  components  are  verified  by  computing  the  same  components,  using 
expressions  derived  by  Wait  [1952,  1961].  The  field  comp'  for  the  IC  are 
calculated  from  explicit  expressions  and  numerical  inte<  ...n  [Wait,  1953;  Inan 
et  al..  1982].  Most  of  the  data  are  presented  in  dimensiomess  form  by  measuring  all 
distances  in  terms  of  the  skin  depth  of  the  sea  water  6,  define  as 

6,  =  (2/u;/io<r,)1/2 

where  jj  is  the  angular  frequency  related  to  frequency  /  by  the  relation  u-—2 irf. 
The  sea-bed  conductivity  is  normalized  to  the  sea-water  conductivity  by  a1  = 
cTf/cr,.  where  a'  is  the  normalized  sea-bed  conductivity.  The  skin  depth  in  the  sea 
bed  is  related  to  the  skin  depth  in  sea  water  by  6f  =  bf/sfip,  and  the  wavelengths 
in  both  media  are  proportional  to  the  skin  depths  because  of  the  relation 

X(*./)  =  27r5(«./) 

In  conducting  media,  the  field  intensities  variy  with  an  exponential  term 
(where  d  is  the  distance  from  the  source)  representing  55  dB/wavelength  attenuation 
[Hansen,  1963],  and  it  is  this  attenuation  that  limits  the  range  of  the  electromagnetic 
signal.  Figure  1.2  is  a  logarithmic  plot  of  skin  depth,  wavelength,  and  exponential 
attenuation  in  sea  water  as  a  function  of  frequency. 

The  numerical  results  obtained  in  this  work  should  be  useful  in  determining 
an  effective  sea-bed  conductivity  in  smooth  regions  of  the  sea  floor.  Stretches  of 
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Figure  1.1.  Geometry  for  a  source  and  receiver  located  on  the  sea  floor. 


the  abyssal  plain  in  the  western  North  Atlantic  have  been  measured  to  be  smooth 
within  2  m  over  distances  of  100  km  [Pickard  and  Emery,  1982).  It  is  also  assumed 
that  the  sea  water  shields  the  sea  floor  from  atmospheric  noise  and,  because  the 
region  of  interest  is  far  from  the  shore,  the  possibility  of  atmospheric  and  power-line 
noise  propagating  through  the  Earth's  crust  can  be  neglected.  This  is  important 
because  atmospheric  noise  is  especially  strong  in  the  ULF/ELF  frequency  range 
(Liebermann,  1962;  Soderberg,  1969],  The  only  other  possible  sources  of  noise  are 
thermal  noise,  motions  of  the  sea  water  (assumed  to  be  negligible),  the  internal 
noise  of  the  receiver,  and,  very  speculatively,  sources  within  the  Earth’s  crust. 


FREQUENCY  (HERTZ) 


Figure  1.2.  Skin  depth,  wavelength,  and  exponential  attenuation  in  sea  water 
a  conductivity  of  4  S/m  (based  on  Figure  A.3  by  Kraichman  (1976]). 


Chapter  II.  Vertical  Magnetic  Dipole 


In  this  chapter  we  consider  the  electromagnetic  fields  produced  by  a  small 
current-carrying  insulated  wire  loop  located  at  the  plane  interface  of  two  semi¬ 
infinite  dissipative  media.  As  illustrated  in  Figure  2.1,  the  loop  is  taken  to  be  at 
the  origin  of  a  cylindrical  coordinate  system  (p,  <(>,  2),  with  its  axis  perpendicular  to 
the  interface  2=0  (vertical  magnetic  dipole).  The  strength  of  the  magnetic  dipole  is 
IdA ,  where  /  is  the  amplitude  of  the  loop  current  Icosut  and  dA  is  the  elementary 
area  of  the  loop. 

A.  Derivation  of  the  Field  Components 

Following  Wait  [1952]  to  derive  expressions  for  the  field  components,  the  mag¬ 
netic  vector  potentials  Ft  and  Ff  are  first  determined  for  the  upper  and  lower 
regions,  respectively.  They  have  only  the  following  2-components: 


where 


F»,z  = 
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iupoIdA  f°° 

J0(\p)\d\ 
u.  +  ti/  uv 
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Here,  Jo(\p)  is  the  Bessel  function  of  the  first  kind  of  zero  order.  The  electric 
and  the  magnetic  fields  are  related  to  these  potentials  by 

E,  =  -  V X  F,  (2.3) 


B,  =  — uo<7iF,-  +  —  W*Fj 
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Figure  2.1.  Vertical  magnetic  dipole  at  the  sea/sea-bed  interface,  (a)  Side  view,  (b) 
top  view,  showing  azimuthal  symmetry. 
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with  E{iP  =  Bi  ^  ==  EitZ  =  0.  The  subscript  »  =  a  or  /  depending  on  whether  the 
fields  are  in  the  upper  or  lower  region;  therefore, 
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(2.6c) 


When  the  observation  point  is  on  the  2-axis  ( p  —  0),  there  is  only  one  nonzero 
field  component  (the  2-component  of  the  magnetic  field)  along  the  2-axis, 


2  JT  Jo  u,  +  Ur 


which  can  also  be  expressed  as 

(27> 

This  expression  can  be  rewritten  in  a  dimensionless  form  by  normalizing  all  dis¬ 
tances  by  the  skin  depth  of  the  sea  water  as 

=  **A.  =  A'  (2.8) 


where 


=  >/x«  +  i2,  +  t'2^ 

with  2/  =  2/6,,  X'  =  6,X,  and  o'  =  <7j/<t,.  Note  that,  when  w=0,  this  magnetic 
field  component  becomes 


B.,z  = 


Ho  Id  A 
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When  the  observation  point  is  on  the  plane  interface  (2=0),  the  field  com¬ 
ponents  in  Eqs.  (2.6)  with  the  exponential  term  being  unity  can  be  expressed  in  a 
dimensionless  form  as 
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where  /  =  p/6,.  Following  Wait(1952],  two  of  the  the  field  components  ( Es ^  and 
have  the  following  explicit  expressions: 

£..♦  =  ..  IT7^,  J3  +  *W>  +  I^V"  -  (3  +  +  l5/>Vw I  (21°l 


=  -  /20/~~2T~5  K9  + 
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(2.11) 


-(9  +  9 qr/p  +  4qf}p2  +  Tfjp3)e  'W] 

The  magnetic-field  component  B,tP  and  the  electric-field  component  E,^  are  zero 
at  oj=0,  and  the  only  nonzero  field  component  B,}Z  at  u»=0  is 


B',t  = 


—Hold A 


B.  Numerical  Results 

The  three  nonzero  field  components  at  the  sea/sea-bed  interface  (2=0)  are 
numerically  evaluated  from  the  integral  expressions  in  Eqs.  (2.6),  based  on  the 
techniques  described  by  Bubenik  [1977].  The  dipole  moment  is  set  equal  to  unity 
(1  Am2)  in  these  computations;  the  field  amplitudes  for  any  other  dipole  moment 
can  be  obtained  quite  simply  from  the  computed  values  by  multiplying  by  the 
appropriate  dipole  moment.  Note  that  the  electric  field  unit  used  in  the  data 
presentation  is  the  microvolt/meter  and  the  magnetic  field  unit  is  the  picotesla. 


The  amplitudes  of  the  field  components  are  plotted  in  a  dimensionless  form  for 
sea-bed  conductivities  ranging  from  1  to  0.001  (1,  0.3,  0.1,  0.03,  0.01,  0.003,  and 
0.001)  times  the  conductivity  of  the  sea  water 

Figure  2.2  shows  the  variations  of  the  magnetic  field  components  Bp  and  Bz, 
total  magnetic  field  Bt ,  and  total  electric  field  Ej{=  Ej,)  as  a  function  of  the 
receiver  distance  normalized  by  the  skin  depth  of  sea  water.  As  indicated,  the  fields 
are  all  presented  in  dimensionless  form.  The  results  are  valid  for  all  frequencies 
in  the  range  over  which  the  assumptions  made  in  the  derivations  are  valid.  The 
actual  field  components  can  be  obtained  by  substituting  a  numerical  value  for  the 
skin  depth  Ss  at  the  frequency  of  interest.  The  observer  distance  ranges  from  0.1 
to  100  sea  water  skin  depths.  Note  that  both  the  vertical  and  horizontal  axes  are 
logarithmic. 

The  ratios  of  the  field  components  produced  on  the  sea  floor  to  the  field  com¬ 
ponents  produced  in  sea  water  of  infinite  extent  (i.e.,  no  sea  floor)  are  plotted  vs 
normalized  distance  in  Figure  2.3.  The  p-component  of  the  magnetic  field  is  zero 
when  there  is  no  sea  floor  present  and,  as  a  result,  no  ratio  curves  are  shown  for 
this  component. 

It  is  also  possible  to  plot  the  field  components  in  another  dimensionless  form 
in  which  the  horizontal  axis  shows  the  frequency  variation  for  a  certain  receiver 
distance.  In  Figure  2.4,  the  field  components  are  parametric  in  terms  of  the  actual 
distance  of  the  receiver,  and  the  plots  reveal  the  frequency  variation  of  these 
parametric  components.  Both  axes  are  drawn  with  logarithmic  scales.  Again,  the 
actual  field  at  any  receiver  distance  can  be  derived  from  the  curves  that  are  shown 
by  substituting  the  value  of  the  receiver  distance. 

The  results  of  the  numerical  integration  were  verified  by  calculating  two  of  the 
field  components  E \  and  Bz ,  using  the  explicit  expressions  in  Eqs.  (2.10)  and  (2.11). 
Some  of  the  low-frequency  portions  of  the  curves  were  also  checked  by  comparing 
their  indicated  field  values  with  those  obtained  from  the  dc  expressions  for  the  field 
components. 


Figure  2.2.  Variations  in  amplitudes  of  the  electric  and  magnetic  fields  produced 
on  the  sea  floor  by  a  vertical  magnetic  dipole  as  a  function  of  distance.  Note 
that  there  are  two  components  of  the  magnetic  field,  one  parallel  and  the  other 
perpendicular  to  the  sea/sea-bed  interface,  and  only  one  component  of  the 
electric  field  parallel  to  the  interface.  The  dipole  moment  is  unity. 
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Figure  2.3.  Variations  in  the  electric  and  the  magnetic  fields  produced  on  the  sea 
floor  for  an  electrically  conducting  sea  bed.  The  curves  show  the  variation  with 
distance  of  the  ratio  of  the  electromagnetic  fields  on  a  sea  floor  to  the  fields  in 
an  infinitely  deep  sea  water. 


C.  Summary 

Based  on  the  numerical  results,  the  presence  of  the  sea  floor  enables  the  fields 
to  propagate  to  longer  distances  because  of  the  assumed  lower  conductivity  of  the 
sea  bed.  To  take  a  particular  illustrative  example,  the  amplitude  of  the  parametric 
electric  field  is  \Epg  J  ~  1.07  at  a  receiver  distance  of  10(5,  on  the  z— 0  plane  in  an 
infinitely  deep  sea.  On  the  sea  floor,  it  reaches  the  same  value  at  a  receiver  distance 
greater  than  156,  for  a  sea-bed  conductivity  of  0.4  S/m,  at  a  distance  greater  than 
226,  for  0.04  S/m,  and  at  a  distance  greater  than  256,  for  0.004  S/m. 

A  new  field  component  Bp  is  produced  as  a  result  of  the  presence  of  the  sea  bed, 
and  it  becomes  larger  than  Bz  at  longer  receiver  distances.  For  example,  at  100  Hz, 
and  for  a  sea-bed  conductivity  of  0.001  times  the  sea-water  conductivity,  \B,  p\  ~ 
6.55  X  10~5  pT  and  —  5.70  X  10-®  pT  at  632  m.  Thus  the  total  magnetic 
field  at  632  m  is  approximately  equal  to  this  new  component.  The  amplitude  of  the 
horizontal  component  of  the  magnetic  field  becomes  \B,tP\  ~  5.70  X  10-6  pT  at 
a  receiver  distance  of  1100  m.  This  new  component  is  also  more  sensitive  to  the 
conductivity  of  the  sea  bed  at  low  frequencies.  For  example,  |£J  p|  =  0.35  X  10®  for 
tr/=0.3,  0.45  X  10®  for  <r'=0.1,  0.49  X  10®  for  <r'==0.01;  the  other  two  components 
are  equal  to  \EP  J  =  0.16  X  108  and  | =  0.10  X  109  for  all  the  above  values 
of  a’ . 

As  can  be  seen  in  Figure  2.4,  the  amplitudes  of  the  fields  at  a  fixed  point  reach 
a  maximum  at  a  certain  optimal  frequency.  Knowledge  of  this  frequency  may  be 
useful  for  two  reasons.  First,  if  the  signal  is  required  to  be  as  strong  as  possible, 
this  frequency  should  be  chosen  (assuming  that  all  frequencies  are  are  equally 
easy  to  generate).  The  disadvantage  is  the  limited  bandwidth,  which  restricts  the 
amount  of  information  transmission  because  the  optimal  frequency  is  very  low. 
Second,  the  effects  of  the  sea  floor  become  significant  around  this  frequency.  In 
Figure  2.4,  the  curves  below  the  optimal  frequency  are  nearly  the  same  for  all  sea¬ 
bed  conductivities;  however,  above  this  frequency,  they  deviate  from  each  other 
depending  on  the  conductivity  of  the  sea  bed. 


Chapter  III.  Horizontal  Electric  Dipole 


In  this  chapter  we  consider  the  electromagnetic  fields  produced  by  an  elementary 
current  source  located  at  the  plane  interface  of  two  semi-infinite  dissipative  media. 
As  illustrated  in  Figure  3.1,  the  source  is  assumed  to  be  located  at  the  origin  of 
a  cylindrical  coordinate  system  (p,  <t>,  2),  with  the  plane  of  separation  of  the  two 
media  corresponding  to  ^=0.  The  source  is  oriented  along  the  x-axis  (».e.,  the  <f>==  0 
axis)  parallel  to  the  interface  (horizontal  electric  dipole)  and  has  an  electric-dipole 
strength  of  /<//,  where  /  is  the  amplitude  of  the  sinusoidal  source  current  I cos(ut) 
and  dl  is  the  elementary  length  of  the  source. 

A.  Derivation  of  the  Field  Components 

For  convenience,  the  Hertz  vector  components  are  used  instead  of  the  field 
components.  Following  Wait  [1961],  the  rectangular  components  of  the  Hertz  vector 
in  the  upper  medium  are 
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and  the  other  terms  are  defined  in  Chapter  II. 

The  electric  and  magnetic  fields  in  the  upper  medium  can  be  obtained  using 
the  two  vector  expressions, 
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where  ne  p,  FJ,  and  /7«,2  are  the  components  in  the  upper  medium  in  cylindrical 
coordinates.  The  divergence  of  the  Hertz  vector  is 
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Based  on  Eqs.  (3.4)  the  field  components  can  be  obtained  in  the  following  form: 
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where  the  argument  of  the  Bessel  functions  Jo  and  J\  are  the  same  as  Eqs.  (3.1). 
The  field  components  can  also  be  written  in  a  dimensionless  form  as  follows: 
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Note  that  the  terms  on  the  left-hand  side  are  independent  of  receiver  distance  so 
that,  after  the  parametric  fields  are  numerically  determined,  they  can  be  useful  for 
all  frequencies  in  the  valid  frequency  range. 

Equations  (3.6)  can  be  simplified  by  considering  two  observer  locations.  The 
first  is  along  the  x-axis  (p  —  0)  where  the  only  nonzero  field  components  are  the 
x-component  of  the  electric  field  E,tX  and  the  j/-component  of  the  magnetic  field 
Bg  y  given  by 
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The  second  observer  is  located  at  the  interface  separating  the  two  media  (2=0), 
where  the  field  components  are  the  same  as  in  Eqs.  (3.6)  but  without  e~u,z  in  all 
the  integrals.  Wait  [1961]  derived  an  explicit  expression  for  the  2-component  of  the 
magnetic  field  B$tZ  in  the  following  form: 
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Three  of  these  field  components  Et,z ,  Bt>p,  and 
the  other  three  simplify  to 

are  zero,  for  the  dc  case,  and 
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B.  Numerical  Results 

The  electric  and  magnetic  field  components  at  the  sea-floor  interface  (z=0) 
are  numerically  evaluated  from  the  integral  expressions  for  two  azimuthal  angles 
4>  =  0°  and  <f>  =  90°,  using  the  techniques  described  by  Bubenik  [1977].  The  field 
components  at  any  arbitrary  <p  can  be  obtained  from  these  results  by  multiplying 
them  by  either  sin  <6  or  cos0.  The  electric  dipole  moment  is  set  equal  to  unity  (1 
Am),  and  the  field  components  for  any  arbitrary  dipole  moment  can  be  determined 
by  multiplying  these  results  by  the  dipole  moment;  the  units  of  the  electric  and 
magnetic  fields  are  in  microvolts/meter  and  picoteslas.  The  results  are  plotted  in 
a  dimensionless  form  for  various  sea-bed  conductivities  ranging  from  1  to  0.001  (1, 
0.3,  0.1,  0.03,  0.01,  0.003,  0.001)  multiplied  by  the  sea-water  conductivity.  The  axes 
are  logarithmic. 

The  amplitudes  of  the  nonzero  field  components  EliP,  Bt>^,  and  E,it  and  the 
total  electric  field  Et  j  are  plotted  vs  distance  in  terms  of  the  skin  depth  of  sea 


water  for  an  azimuthal  angle  of  <j>  =  0°,  as  illustrated  in  Figure  3.2.  Figure  3.3 
plots  the  variations  of  the  amplitudes  of  the  nonzero  field  components  BtiP,  E 
and  B,'Z  and  the  total  magnetic  field  B,  j  when  <j>  =  90°.  All  the  results  can  be 
converted  to  the  real  field  values  for  any  frequency  by  substituting  the  numerical 
value  of  the  skin  depth  6 ,  at  that  frequency  and  the  conductivity  a,  of  the  sea 
water.  The  receiver  distance  varies  from  0.1  to  100  sea  water  skin  depths. 

The  numerical  results  for  the  field  components  in  the  alternate  dimensionless 
form  can  be  obtained  by  holding  the  distance  of  the  observer  constant  and  varying 
the  frequency,  as  demonstrated  in  Figures  3.4  and  3.5,  where  the  increase  in  the 
horizontal  axis  shows  the  rise  in  frequency.  The  actual  field  values  for  any  arbitrary 
receiver  distance  p  can  be  calculated  from  these  curves  by  substituting  the  numerical 
value  of  p  in  the  parametric  expressions. 

The  2-component  of  the  magnetic  field  BtfZ  was  also  calculated  from  the  explicit 
expression  in  Eq.  (3.10),  and  the  results  were  used  to  verify  those  obtained  by  the 
numerical  integration  of  B,iZ.  The  dc  expressions  for  some  of  the  field  components 
were  also  used  to  check  the  lower  frequency  portion  of  the  curves. 

C.  Summary 

The  ranges  of  the  electromagnetic  signals  tend  to  increase  as  a  result  of  the 
lower  conductivity  of  the  sea  bed,  as  indicated  by  all  the  numerical  results.  For 
example,  when  the  sea-bed  conductivity  is  0.01  times  the  sea-water  conductivity, 
the  increase  is  approximately  a  factor  of  105  at  distances  of  10  to  20  skin  depths. 

The  new  field  components  produced  by  the  existence  of  the  sea  floor  are  the  p- 
and  (^components  of  the  magnetic  field  BtiP  and  Bti^  and  the  2-component  of  the 
electric  field  Et<z.  The  new  magnetic  field  components  B$tP  and  B, ^  become  larger 
than  the  2-component  of  the  magnetic  field  Bt>z  at  longer  distances,  which  results 
in  a  greater  total  magnetic  field.  The  new  component  of  the  electric  field  does  not 
have  this  property;  however,  because  of  its  sensitivity  to  sea-bed  conductivity  even 
at  low  frequencies,  it  may  have  useful  applications  in  sea-bed  prospecting. 

Figures  3.4  and  3.5  show  the  variations  of  the  fields  with  frequency.  As  can  be 
seen,  all  of  the  field  components  are  sensitive  to  the  conductivity  of  the  sea  bed 
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Figure  3.2.  Variations  in  the  amplitudes  of  the  electric  and  magnetic  fields  produced 
on  the  sea  floor  by  a  horizontal  electric  dipole  at  an  azimuthal  angle  of  0°  as  a 
function  of  distance.  Note  that  there  are  two  components  of  the  electric  field,  one 
parallel  and  the  other  perpendicular  to  the  sea  floor,  and  only  one  component 
of  the  magnetic  field,  directed  parallel  to  the  floor.  The  dipole  moment  is  unity. 
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Figure  3.3.  Variations  in  tlie  amplitudes  of  the  electric  and  magnetic  fields  produced 
on  the  sea  floor  by  a  horizontal  electric  dipole  at  an  azimuthal  angle  of  90°  as 
a  function  of  distance.  Note  that  there  are  two  components  of  the  magnetic 
field,  one  parallel  and  the  other  perpendicular  to  the  sea  floor,  and  only  one 
component  of  the  electric  field,  directed  parallel  to  the  floor.  The  dipole  moment 


at  low  frequencies  except  for  the  vertical  component  of  the  magnetic  field.  The 
vertical  component  of  the  electric  field  is  interesting  because  it  is  highly  sensitive  to 
sea-bed  conductivity  at  low  frequencies  and  has  an  optimal  frequency  that  is  also 
sensitive  to  the  same  conductivity.  Unlike  the  magnetic  field  components  B,%f  and 
B, produced  at  the  sea  floor,  this  component  is  not  large  at  long  distances  and, 
as  a  result,  it  does  not  contribute  to  the  total  electric  field. 


Chapter  IV.  Infinite  Cable 


In  this  chapter  we  consider  the  electromagnetic  fields  produced  by  a  straight 
current-carrying  insulated  cable  of  infinite  length  lying  at  the  plane  interface  of 
two  conducting  media — in  our  case  the  interface  between  the  sea  and  the  sea  bed. 
We  use  the  term  “infinite”  in  a  relative  sense:  if  the  length  of  the  cable  is  much 
greater  than  all  the  other  relevant  linear  dimensions,  it  is  known  as  an  “infinite” 
cable.  It  is  assumed  to  carry  an  equally  distributed  alternating  current  Icosut  at  a 
given  instant  of  time,  and  this  assumption  is  justified  at  sufficiently  low  frequencies 
[Wait,  1952a;  Sunde,  1968].  The  cable  is  oriented  along  the  x-axis,  and  the  plane 
interface  of  the  two  media  is  2=0.  Figure  4.1  shows  the  geometry  in  the  y,  z  plane; 
everything  is  invariant  in  the  x-direction  and  the  direction  of  the  current  is  taken 
to  be  out  of  the  page  at  /=0. 


A.  Derivation  of  the  Field  Components 

The  Hertz  vector  components  in  the  upper  medium  for  an  HED  lying  at  x=0, 
oriented  along  the  x-axis,  and  located  on  the  x-axis  at  x  =  l  are  the  same  as  in  Eq. 
(3.1)  with  x  replaced  by  x  —  l.  For  an  infinite  cable,  the  Hertz  vector  components 
can.be  derived  by  integrating  the  components  of  the  Hertz  vector  for  an  HED  from 
/  =  —oo  to  /  =  oo,  which  results  in 
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and  all  the  other  quantities  were  defined  in  Chapter  HI.  The  partial  derivative  d/dx 
can  be  replaced  by  —d/dl  to  yield 


Because  the  divergence  of  the  Hertz  vector  is  also  zero,  the  only  nonzero  vector 
component  is  the  .r-component.  By  using  the  integral  [Sunde,  1968] 


the  expression  for  this  component  can  be  simplified  to 
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The  electric  and  magnetic  field  components  can  be  derived  from  the  Hertz  vector 
via  Eqs.  (3.2)  and  (3.3)  as 
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where  Btx  =  ESiy  =  =  0. 

As  in  Chapters  II  and  HI,  two  receiver  positions  are  considered  which,  to  some 
extent,  will  simplify  the  above  expressions.  When  the  receiver  is  located  along  the 
line  perpendicular  to  the  interface  and  passes  through  the  infinite  cable  (y=0), 
B,  .=  0  and  the  other  two  field  components  become 
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These  expressions  can  be  further  simplified  by  multiplying  the  numerator  and 
denominator  of  the  integrands  by  u,  —  ur  and  following  an  approach  similar  to 


that  used  by  Wait  and  Spies  [1971];  then, 
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Note  that  A'ohgz)  is  the  modified  Bessel  function  of  the  second  kind  of  zero  order 
and  with  complex  argument.  In  the  dc  case,  the  above  two  components  can  be 
written 
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When  the  receiver  is  located  at  the  interface  2=0,  the  field  components  can  be 
written  in  dimensionless  form  as 
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(4.6c) 


where 

u't  =  \J\n  +  *2,  «y.  =  i/\'2  +  1 2^ 

and  2;  =  2/6,,  y*  =  y/6,,  X'  =  X6,,  and  <7*  =  Of/o,.  Wait  [1953,  1962]  derived 
explicit  expressions  for  the  horizontal  electric  field  component  Etx  and  the  vertical 
magnetic  field  component  Bt  2  as 


EtI 


xu  Hoi 

x(Tf-i2,)y2 


[i»yKih,y)  -  ifyKi(i/y)\ 


(4.7) 


33 


B,,z  = 
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where,  again,  K o  and  K\  are  the  modified  Bessel  functions  of  the  second  kind 
of  order  zero  and  one  with  complex  arguments.  The  horizontal  component  of  the 
magnetic  field  B,  y  does  not  have  an  explicit  expression.  When  cj=0,  both  E,iX  and 
Bg  y  become  zero  and  the  vertical  magnetic  field  component  simplifies  to 
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B.  Numerical  Results 

This  section  presents  the  numerical  data  for  the  electric  and  magnetic  field 
components  at  the  interface  of  the  two  conducting  media.  Two  of  these  components, 
Et  l  and  Btz,  have  explicit  expressions  [Eqs.  (4.7)  and  (4.8)]  that  can  be  expressed  in 
terms  of  Kelvin  functions  (Young  and  Kirk,  1964]  and  their  derivatives.  After  some 
algebraic  manipulations,  their  amplitudes  can  be  written  in  dimensionless  form  as 

\EU  =  ’^J1\E.A  =  -  Re*)2  +id>  -  H-»> 

and 

\BP,,\  =  ~~ \B.A  =  (1 .  ff2)Q2  K".(«i  -  -  2(<*i  -  **))* 

(410) 
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where 


a\  =ker0at, 
=kei0at, 
ci  =ker'0a„ 
d\  =kei'0a„ 


a2  —  ker0a  f 
b2  =  kei^otf 
c2  =  ker oa  / 
d2  =  kei'Qa  f 
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»  . 


and 

ae  =  v/2  j/,  a;  =  Rag}  R2  =  o'  =  <ry/<7, 

where  ae  and  ay  are  defined  as  the  induction  numbers  in  the  sea  water  and  sea 
bed  respectively  [Coggon  and  Morrison,  1970].  The  above  expressions  for  Ee,x  and 
B,,x  can  be  evaluated  numerically  by  using  tables  of  Kelvin  functions  and  their 
derivatives  [Lowell,  1959]. 

The  other  nonzero  component  at  the  interface  is  the  horizontal  magnetic  field 
component  Be  y,  and  it  can  be  computed  only  by  numerical  integration.  The  ex¬ 
pression  for  it  [Eq.  (4.6b)]  can  be  divided  into  two  parts, 
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Two  factors  must  be  determined  for  numerical  integration.  The  first  is  \'max.  It  is 
assumed  that  the  effective  conductivity  of  the  sea  floor  is  smaller  than  that  of  sea 
water  (o'  <  1).  The  value  for  X'max  is  chosen  such  that  X^0I  >>  2  and,  if  this 
constraint  is  satisfied,  both  u't  and  u'y  can  be  approximated  as 
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for  X'  >  \'max.  Based  on  this  approximation,  the  second  integral  above  can  be 
replaced  by  an  explicit  term  yielding 


|4U) 

Similarly  the  integral  expression  for  the  vertical  component  of  the  magnetic  field 
[Eq.  (4.6c)]  can  be  written  as 
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(4.12) 


The  above  integrands  can  be  separated  into  their  real  and  imaginary  parts  and  thus 
numerically  integrated  as  real  integrals.  The  data  for  the  field  components  can  then 
be  obtained  by  recombining  these  terms  according  to  Eqs.  (4.11)  and  (4.12). 
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The  second  factor  that  must  be  determined  is  the  method  of  integration.  Both 
integrands  are  well-behaved  functions  except  for  the  cosine  or  sine  term  that  may 
oscillate  rapidly  at  greater  distances  than  the  skin  depth  of  the  seawater  [Hermance 
and  Peltier,  1970].  Two  techniques  of  numerical  integration,  Weddle’s  rule  [Computation 
Laboratory  of  Harvard  University,  1949]  and  Filon’s  method  [Tranter,  1956],  were 
used.  Provided  the  integration  interval  is  made  small  enough,  both  techniques 
produce  sufficiently  accurate  results. 

The  computations  can  be  verified  in  several  ways.  The  first  one  is  to  compare 
the  results  of  the  two  integration  techniques.  The  second  is  to  compare  the  results 
obtained  for  BtiZ  through  numerical  integration  using  Eq.  (4.12)  to  the  results 
obtained  for  the  same  component  using  Eq.  (4.10)  and  tables  of  Kelvin  functions. 

The  third  is  to  check  the  lower  and  higher  frequency  regions  by  approximate 
expressions  for  these  regions. 

The  left  hand  panels  in  Figures  4.2,  4.3,  4.4,  and  the  single  panel  in  Figure 
4.5  show  the  variations  in  the  parametric  amplitudes  of  the  three  nonzero  field 
components  Ept  x,  B?y,  and  BPe  z  and  the  total  magnetic  field  Bpe  T  with  the  induction 
number  of  sea  water  a*.  Note  that  the  parametric  terms  ncrt6“/I  and  \f2x6tlnQl 
appearing  on  the  left  hand  side  depend  on  frequency  /  and  not  on  receiver  distance 
y  [see  Eqs.  (4.6)].  These  curves  can  be  used  for  any  frequency  in  the  valid  frequency 
range.  The  horizontal  axis  indicates  the  perpendicular  distance  of  the  receiver  from 
the  axis  of  the  cable  along  the  sea  floor  in  terms  of  skin  depths  of  sea  water.  Each 
curve  corresponds  to  a  different  value  of  R  (R  =  y 4ry/<r,);  R  takes  the  values  of 
0.01,  0.03,  0.1,  0.3,  and  1  except  for  the  BPty  component  which  is  zero  when  /?=1 
(no  sea  floor).  Both  axes  are  plotted  on  a  logarithmic  scale,  and  the  units  are  in 
volts/meter  and  teslas. 

The  right  hand  panels  in  Figures  4.2  and  4.4  show  the  variations  of  the  amplitude 
ratio  of  E,yX  and  B,  z  produced  on  the  sea  floor  to  those  produced  in  a  sea  of  infinite 
depth  as  a  function  of  distance.  The  horizontal  magnetic  field  component  BtiS/  is  zero 
in  the  absence  of  the  sea  floor.  The  second  panel  in  Figure  4.3  plots  the  variation  of 
the  ratio  of  the  amplitude  of  to  the  amplitude  of  BtiZ  as  a  function  of  distance. 


Figure  4.2.  Variation  with  horizontal  distance  of  (1)  the  amplitude  of  the  electric 
field  produced  on  the  sea  floor  by  an  infinite  cable  (left  panel)  and  (2)  of  the  ratio 
curves  illustrating  the  changes  produced  in  the  electric  field  by  the  presence  of 
an  electrically  conducting  sea  bed  (right).  Each  curve  corresponds  to  a  different 
sea-bed  conductivity.  The  curves  in  the  right  hand  panel  show  the  ratio  of  the 
electric  field  at  the  sea  floor  to  the  electric  field  in  an  infinitely  deep  sea. 
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Figure  4.3.  Variation  with  horizontal  distance  of  (1)  the  amplitude  of  the  horizontal 
component  of  the  magnetic  field  produced  on  the  sea  floor  by  an  infinite  cable 
(left  panel)  and  (2)  of  the  ratio  of  the  horizontal  component  to  the  vertical 
component  of  the  magnetic  field,  both  components  being  measured  at  the  same 
location  on  the  sea  floor  (right).  Each  curve  is  a  function  of  distance  and 
corresponds  to  a  different  sea-bed  1  ductivity. 


Figure  4.4.  Variation  with  horizontal  distance  of  (1)  the  amplitude  of  the  verti¬ 
cal  component  of  the  magnetic  produced  on  the  sea  floor  by  an  infinite  cable 
(left  panel)  and  (2)  of  the  ratio  curves  illustrating  the  changes  in  the  vertical 
component  of  the  magnetic  field  produced  by  the  presence  of  an  electrically 
conducting  sea  bed  (right).  Each  curve  corresponds  to  a  different  sea-bed  con¬ 
ductivity.  The  curves  in  the  right  hand  panel  show  the  ratio  of  the  vertical 
component  of  the  magnetic  field  produced  on  the  sea  floor  to  the  same  com¬ 
ponent  in  an  infinitely  deep  sea. 
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Figure  4.6.  Variation  with  horizontal  distance  of  the  amplitudes  of  the  horizontal 
and  vertical  components  of  the  magnetic  field  produced  on  the  sea  floor  by  an 
infinite  cable  carrying  an  alternating  current  with  an  amplitude  of  1000  A  and 
a  frequency  of  1  Hz.  Each  curve  corresponds  to  a  different  sea-bed  conductivity. 
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Figure  4.7.  Variation  with  horizontal  distance  of  the  amplitudes  of  the  total  electric 
and  magnetic  fields  produced  on  the  sea  floor  by  an  infinite  cable  carrying  an 
alternating  current  with  an  amplitude  of  1000  A  and  a  frequency  of  1  Hz.  Each 
curve  corresponds  to  a  different  sea-bed  conductivity. 


In  Figures  4.6  and  4.7,  the  real  amplitudes  of  the  field  components  EgtX,  B,iy , 
Be,z  and  the  total  magnetic  field  Bsj  are  plotted  vs  distance  for  an  alternating  cur¬ 
rent  source  with  an  amplitude  of  1000  A  and  a  frequency  of  1  Hz  for  the  above  sea¬ 
floor  conductivities.  The  units  of  electric  and  magnetic  fields  are  millivolts/meter 
and  picotesla.  and  distance  is  in  meters. 

The  above  results  can  be  plotted  in  another  parametric  form  in  which  the 
amplitude  of  the  field  components  is  written  as 
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Here,  the  terms  on  the  left  hand  side  are  not  functions  of  frequency  but,  instead,  are 
functions  of  receiver  distance  y.  Figures  4.8  and  4.9  plot  the  variations  of  the  field 
components  in  Eqs.  (4.13),  and  the  horizontal  axis  shows  the  variation  in  frequency. 
Note  that  each  curve  now  corresponds  to  a  different  sea-floor  conductivity;  again, 
there  is  no  horizontal  magnetic  field  component  for  /2=L 


C.  Summary 

All  three  field  components  produced  at  the  interface  can  propagate  farther  when 
a  sea  floor  is  present  than  can  those  produced  in  the  absence  of  a  sea  floor  because 
of  the  assumed  lower  conductivity  of  the  sea  bed.  The  increase  in  range  becomes 
significant  at  greater  distances.  For  example,  when  /=1  Hz,  the  amplitude  of  the 
horizontal  component  of  the  electric  field  E',*  becomes  1  nV/m  at  a  distance  of 
y  ~  3  km  for  n‘—\  (infinitely  deep  sea),  y  ~  18  km  for  <r>= 10-2,  and  y  ~  90  km 
for  <r#=10“4.  At  90  km,  the  source  cable  must  be  much  longer  for  it  to  be  an  infinite 
cable;  for  example,  if  it  is  ten  times  longer  (900  km),  a  uniform  current  distribution 
is  still  a  good  assumption  at  1  Hz  [Inan  et  al.,  1983]. 

The  new  magnetic  field  component  Bg>y  is  parallel  to  the  interface  and  perpen¬ 
dicular  to  the  axis  of  the  source.  It  is  negligible  when  compared  to  the  vertical 
component  BttZ  near  the  cable;  however,  it  becomes  comparable  to  and  larger  than 
Bt  s  at  greater  receiver  distances  (Figure  4.3).  For  example,  if  7=1000  A,  /=  1 
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Figure  4.9  Variations  in  the  amplitudes  of  the  vertical  component  of  the  magnetic 
field  and  the  total  magnetic  field  produced  on  the  sea  floor  by  an  infinite 
cable  as  a  function  of  frequency.  Each  curve  corresponds  to  a  different  sea-bed 
I  conductivity. 
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Hz.  and  a'  =  0.01,  the  amplitudes  of  the  horizontal  and  vertical  magnetic  field 
components  would  be  \Bg:y\  —  0.70  x  108  pT  and  |B#1.|  cz  0.20  X  108  pT  at  10  m, 
0.57  X  10®  and  0.19  x  107  pT  at  100  m,  and  0.68  X  105  and  0.30  X  105  pT  at  1 
km.  The  new  component  is  also  more  sensitive  to  the  conductivity  of  the  sea  bed 
than  are  the  other  two  field  components  at  lower  frequencies. 

The  horizontal  components  of  the  magnetic  and  electric  fields  are  zero  when 
u.'— —0.  They  also  have  optimal  frequencies,  as  observed  in  Figure  4.8.  For  example, 
at  a  distance  of  1  km,  the  electric  field  component  has  an  optimal  frequency  of 
approximately  0.15  Hz  in  an  infinitely  deep  sea  [Inan  et  al.,  1983]. 


Chapter  V.  Conclusions  and  Recommendations 


A.  Conclusions 

This  investigation  has  concentrated  on  the  electromagnetic  fields  produced 
along  the  sea/sea-bed  interface  (considered  to  be  a  plane  boundary  between  two 
semi-infinite  conducting  media)  by  several  different  harmonic  sources:  a  vertical 
magnetic  dipole,  a  horizontal  electric  dipole,  and  an  infinite  cable.  New  expressions 
and  numerical  results  have  been  obtained  for  the  electromagnetic  fields  produced  by 
these  sources.  The  displacement  current  terms  in  both  media  have  been  neglected, 
which  is  valid  for  frequencies  of  less  than  100  kHz.  The  results  have  been  expressed 
in  two  dimensionless  forms  that  make  them  useful  for  all  frequencies  and  receiver 
distances  in  the  above  frequency  range.  These  new  results  can  be  summarized  as 
follows. 

1)  Vertical  magnetic  dipole:  the  three  nonzero  field  components  at  the  sea/sea¬ 
bed  interface  are  the  horizontal  magnetic  field  component  in  the  p-direction,  the 
horizontal  electric  field  component  in  the  ^-direction,  and  the  vertical  magnetic 
field  component  in  the  ^-direction.  The  first  is  a  new  component  (compared  with 
what  would  be  produced  by  the  same  source  in  a  sea  of  infinite  extent  under 
otherwise  identical  conditions)  that  can  be  useful  because  it  (1)  becomes  larger 
than  the  vertical  magnetic  field  component  at  greater  receiver  distances,  (2)  has 
an  optimal  frequency  at  which  the  amplitude  of  the  field  has  a  maximum,  and  (3) 
is  more  sensitive  to  the  conductivity  of  the  sea  bed  than  the  other  components 
at  low  frequencies.  The  horizontal  electric  dipole  also  has  an  optimal  frequency. 
Both  horizontal  components  are  zero  in  the  dc  case  and  become  more  sensitive  to 
the  conductivity  of  the  sea  bed  above  the  optimal  frequency  than  does  the  vertical 
component  of  the  magnetic  field,  which  is  nonzero  in  the  dc  case. 

2)  Horizontal  electric  dipole:  all  six  field  components  are  nonzero  on  the  sea 
floor  interface.  The  two  azimuthal  angles  considered  are  0  =  0°  and  4>  =  90°. 
At  0  =  0°,  the  three  nonzero  field  components  are  the  horizontal  component  of 
the  electric  field  along  the  axis  of  the  source,  the  horizontal  component  of  the 
magnetic  field  perpendicular  to  the  axis  of  the  source,  and  the  vertical  component 
of  the  electric  field  perpendicular  to  the  sea  floor.  Both  the  horizontal  magnetic 


field  component  and  the  vertical  electric  field  component  are  the  result  of  the 
sea/sea-bed  interface.  The  horizontal  and  especially  the  vertical  component  of 
the  electric  field  are  sensitive  to  the  conductivity  of  the  sea  bed  even  at  very 
low  frequencies.  Another  interesting  feature  is  that  it  has  an  optimal  frequency 
which,  unlike  other  components,  is  also  very  sensitive  to  the  conductivity  of  the  sea 
bed.  At  0  =  90°.  the  three  nonzero  field  components  aie  the  horizontal  magnetic 
field  component  perpendicular  to  the  axis  of  the  source,  the  horizontal  electric 
field  component  parallel  to  the  axis  of  the  source,  and  the  vertical  magnetic  field 
component  perpendicular  to  the  sea  floor.  Only  the  horizontal  component  of  the 
magnetic  field  is  produced  as  a  result  of  the  sea/sea-bed  interface;  it  becomes 
larger  than  the  vertical  component  at  greater  receiver  distances,  is  sensitive  to  the 
conductivity  of  the  sea  bed  at  low  frequencies,  and  is  zero  at  ui— 0.  The  horizontal 
component  of  the  electric  field  at  <f>  —  90°  is  also  sensitive  to  the  conductivity  of 
the  sea  bed  at  low  frequencies. 

3)  Infinite  cable:  the  three  nonzero  field  components  at  the  interface  are  the 
horizontal  electric  field  component  aligned  with  the  axis  of  the  source,  the  horizontal 
magnetic  field  component  perpendicular  to  the  axis  of  the  source,  and  the  vertical 
magnetic  field  component  perpendicular  to  the  sea  floor.  The  horizontal  component 
of  the  magnetic  field  becomes  larger  than  the  vertical  component  at  large  receiver 
distances,  has  an  optimal  frequency,  it  is  sensitive  to  the  conductivity  of  the  sea  bed 
at  low  frequencies,  and  is  zero  at  w=0  The  horizontal  component  of  the  electric 
field  is  also  zero  at  u;=0  and  has  an  optimal  frequency.  The  vertical  component  of 
the  magnetic  field  is  nonzero  even  when  u>=0. 

The  electromagnetic  signals  produced  along  the  sea  floor  can  propagate  farther 
because  of  the  assumed  lower  conductivity  of  the  sea  bed.  Also,  new  field  com¬ 
ponents  are  produced  as  a  result  of  the  sea/sea-bed  interface  and,  because  they 
are  comparable  to  and  even  larger  than  the  components  existing  in  the  infinitely 
deep  sea,  they  produce  larger  total  fields.  For  example,  for  an  infinite  cable  carrying 
alternating  current  with  an  amplitude  of  1000  A  and  a  frequency  of  l  Hz,  for  a  sea¬ 
bed  conductivity  of  1  S/m,  the  amplitude  of  the  vertical  component  of  the  magnetic 
field  is  approximately  0.47  pT  at  a  receiver  distance  of  5  km.  The  amplitude  of  the 
horizontal  component  of  the  magnetic  field  at  5  km  is  approximately  0.75  pT,  which 
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increases  the  amplitude  of  the  total  magnetic  field  to  roughly  0.88  pT.  The  new  field 
components  are  also  sensitive  to  the  conductivity  of  the  sea  bed;  for  example,  in 
the  case  an  infinite  cable,  the  parametric  amplitude  of  the  horzontal  the  magnetic 
dipole  is  0.37,  0.50,  and  0.61  for  o'  =  0.25,  0.09,  and  0.01  at  a  sea-water  induction 
number  of  ae  —  0.1,  whereas  the  parametric  amplitudes  of  the  vertical  magnetic 
dipole  at  the  same  sea-bed  conductivities  and  sea-water  induction  number  are  10. 
Even  this  component  becomes  less  sensitive  to  sea-bed  conductivity,  however,  as 
the  conductivity  is  reduced.  The  optimal  frequency  may  have  some  applications 
for  communicating  at  short  ranges  and  for  sea-bed  prospecting.  For  example,  for  a 
vertical  magnetic  dipole,  the  optimal  frequency  occurs  when  the  receiver  distance  is 
3  to  4  times  the  skin  depth  of  sea  water  and  the  maximum  value  of  the  amplitudes 
of  the  horizontal  magnetic  and  electric  fields  are  sensitive  to  the  conductivity  of 
the  sea  bed  at  low  sea-bed  conductivities. 

B.  Applications 

An  active  experiment  has  been  conducted  on  the  sea  floor  in  the  vicinity  of 
the  East  Pacific  Rise  to  investigate  the  conductivity  of  the  sea  bed  by  analyzing 
the  measured  data  [Young  and  Cox,  1981].  The  numerical  results  that  have  been 
obtained  in  Chapters  II,  III,  and  IV  could  be  applied  to  investigate  the  effective 
conductivity  of  the  sea  bed  (assumed  to  be  homogeneous).  Two  methods  can  be 
used.  One  is  to  keep  the  frequency  of  the  source  constant  and  to  vary  the  observation 
distance  along  the  interface;  however,  most  of  the  field  components  become  sensitive 
to  the  effective  conductivity  of  the  sea  bed  at  distances  greater  than  a  skin  depth  of 
sea  water.  The  other  is  to  hold  the  location  of  the  observer  constant  at  the  interface 
and  to  vary  the  frequency  of  the  source;  in  this  case  the  optimal  frequency  may 
yield  useful  information. 

Another  application  is  the  use  of  electromagnetic  signals  for  undersea  com¬ 
munication.  The  range  of  communication  with  submersibles  in  the  deep  parts  of 
the  ocean  is  limited  because  of  the  high  conductivity  of  the  sea  water.  The  sea 
bed  has  an  effective  conductivity  smaller  than  that  of  sea  water  and,  as  a  result,  it 
may  provide  a  new  path  for  the  signals  along  which  there  will  be  less  attenuation 
compared  to  a  signal  propagating  directly  through  sea  water;  longer  ranges  can  be 


achieved  before  the  signal  becomes  weak  and  cannot  be  detected  because  of  external 
and  internal  receiver  noise.  The  sea  bed  is  also  a  conducting  medium  and  exponen¬ 
tial  attenuation  is  still  significantly  large  at  high  frequencies;  therefore,  operating 
frequencies  must  be  low  enough  to  reach  long  ranges.  This  limits  the  bandwidth, 
and  information  cannot  be  transferred  at  a  high  data  rate. 

New  field  components  at  the  sea/sea-bed  interface  will  increase  the  ranges  of 
communication.  For  an  infinite  cable  carrying  an  alternating  current  of  1000  A  at 
1  Hz,  and  assuming  the  minimum  measurable  magnetic  field  to  be  0.1  pT,  the  field 
on  the  sea  floor  can  be  detected  over  a  range  of  approximately  4  km  before  its 
amplitude  drops  below  0.1  pT.  The  amplitudes  of  the  new  horizontal  component 
and  the  vertical  component  of  the  magnetic  field  along  the  interface  are  0.24  and 
0.09  pT  at  a  receiver  distance  of  8  km  for  an  effective  sea  bed-conductivity  of  0.4 
S/m,  thus  producing  a  total  magnetic  field  of  0.26  pT.  If  the  sea-bed  conductivity 
is  0.04  S/m,  then  the  two  field  components  are  0.27  and  0.03  pT  at  a  receiver 
distance  of  20  km,  and  the  total  magnetic  field  becomes  0.27  pT  as  a  result  of  this 
new  component.  These  two  components  are  about  0.39  X  103  and  0.22  X 102  pT  at 
20  km  for  a  sea-bed  conductivity  of  0.004  S/m.  These  examples  indicate  that  the 
range  of  electromagnetic  fields  along  the  sea  floor  increases  significantly  before  their 
amplitudes  drop  below  the  minimum  measurable  field  value. 

Arrays  of  long  cables  located  on  the  sea  floor  can  achieve  longer  communication 
ranges  [Inan  et  al.,  1982].  The  phases  of  the  currents  in  each  cable  can  be  adjusted 
to  produce  a  maximum  field  amplitude  at  the  receiver,  which  makes  communication 
possible  over  a  much  larger  area.  The  separation  distance  of  the  cables  depends  on 
the  amplitudes  of  the  current,  frequency,  conductivity  of  the  sea  bed,  sensitivity 
of  the  receiver,  and  background  noise.  With  a  series  array  of  cables,  each  carrying 
an  alternating  current  of  amplitude  1000  A  and  frequency  1  Hz  and  separated 
from  each  other  by  30  skin  depths,  the  magnetic  field  produced  at  a  point  midway 
between  the  two  cables  would  be  approximately  4  X  102/£,  or  1.6  pT  in  an  infinitely 
deep  sea;  on  a  sea  bed  with  <7y=0.04  S/m  this  field  would  be  4.6  X  103  pT,  which 
indicates  that  the  cables  could  be  moved  apart  from  a  spacing  of  30£,  to  140<5»  while 
still  yielding  the  same  field  at  the  center  point.  Large  amounts  of  power  will  be  lost 
due  to  resistive  heating  in  both  cases,  however,  as  the  source  current  flows  in  both 
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the  long  cables  and  surrounding  conducting  media.  The  advantage  of  this  system 
is  a  less  noisy  environment  in  which  the  electromagnetic  signals  propagate,  since 
there  is  significant  shielding  of  atmospheric  noise  by  the  large  bulk  of  sea  water 
[Mott  and  Biggs,  1963]. 

C.  Recommendations 

Throughout  this  investigation,  it  has  been  assumed  that  the  sea  bed  is  a  homo¬ 
geneous  medium;  however,  it  actually  contains  layers  of  sediment  with  different 
properties.  This  work  could  be  extended  to  two  or  more  layers.  For  example,  in  the 
two-layer  bed  in  Figure  5.1,  the  first  layer  has  a  conductivity  of  a/i  and  extends 
from  2=0  to  z=—h  and  the  second  has  a  conductivity  of  <jj%  and  extends  from 
z=—h  to  2=— oo.  Both  the  source  and  receiver  are  located  at  the  sea/sea-bed 
interface.  The  conductivity  of  the  upper  layer  is  larger  than  that  of  the  lower  layer. 
For  receiver  distances  smaller  than  the  depth  of  the  first  layer,  the  effects  of  the 
low-conductivity  layer  should  not  be  very  significant.  As  receiver  distances  become 
larger  than  the  depth  of  the  first  layer,  however,  the  signal  that  follows  the  down- 
under-up  mode  will  be  much  stronger  than  the  signal  that  directly  propagates  from 
the  source  to  the  receiver  using  the  sea/sea-bed  boundary,  and  the  lower  layer  will 
be  important.  This  means  that,  depending  on  the  location  of  the  observer,  the  signal 
at  the  receiver  will  carry  information  about  the  properties  of  different  layers  of  the 
sea  bed. 

The  length  of  the  straight  current-carrying  cable  source  was  assumed  to  be 
infinite.  This  is  a  good  assumption  for  a  receiver  located  in  the  vicinity  of  the  cable 
(i.e.,  the  perpendicular  distance  to  the  receiver  is  much  smaller  than  the  length  of 
the  cable)  and  not  close  to  the  ends  of  the  cable.  A  uniform  current  distribution 
along  the  cable  was  also  assumed.  Although  this  is  not  true  in  an  infinite  cable, 
the  fields  produced  at  the  receiver  will  be  the  result  of  the  middle  portion  of  the 
cable  because  the  contribution  from  the  portions  toward  the  ends  (not  necessarily 
carrying  the  same  current)  will  be  negligible.  As  the  receiver  distance  increases, 
however,  a  longer  portion  of  the  cable  must  be  considered  and  the  uniform-current 
assumption  may  no  longer  be  valid.  One  extension  of  this  work  would  be  to  calculate 
the  fields  produced  on  the  sea  floor  by  a  straight  current-carrying  insulated  cable 
of  finite  length. 
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Figure  5.1.  Geometry  for  a  sea  bed  containing  two  conducting  layers.. 


In  practice,  it  may  not  be  possible  to  locate  the  receiver  precisely  at  the  interface. 
Another  extension  would  be  to  place  the  receiver  some  distance  above  the  sea/sea¬ 
bed  interface  and  then  calculate  the  fields  for  this  configuration.  The  integrals  for 
the  field  expressions  will  then  require  an  additional  exponential  term  e~u,z  because 
of  the  attenuation  of  the  signal  as  it  propagates  the  vertical  distance  between  the 
source  and  receiver  in  the  sea  water.  As  a  result  of  this  attenuation,  the  signals  at 
some  vertical  distance  away  from  the  source  will  be  smaller  than  on  the  sea  floor. 
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H.  Ko  I 

Johns  Hopkins  Road 
Laurel,  MD  20810 

La  Jolla  Institute 
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La  Jolla,  C'A  92407 
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P.O.  Box  398 
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University  of  Texas,  Austin 
Geomagnetics  and  Electrical 
Geoscience  Laboratory 
ATTN:  F.X.  Bostick,  Jr.  1 
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Commander 
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University  of  California 
Scripps  Institute  of  Oceanography 
ATTN:  C.S.  Cox  1 

La  Jolla,  CA  92093 

Lockheed  Palo  Alto  Research 
Laboratory 

ATTN:  J.B.  Cladis  1 

W.L.  Imhof  1 

J.B.  Reagan  1 

M.  Walt  1 

3251  Hanover  Street,  Bldg  255 
Palo  Alto,  CA  94304 


Chief 
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Commander 

Air  Force  Systems  Command  2 

Andrews  AFB,  MD  20331 

Commander 

Air  Force  Geophysics  Laboratory 
ATTN:  Technical  Library  1 

Hanscom  AFB,  MA  01731 


m 


END 


2-85 


DTIC 


